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We analyze the simplest possible realization of the curvaton scenario, where a nearly scale- 
invariant spectrum of adiabatic perturbations is generated by conversion of an isocurvature per- 
turbation generated during inflation, rather than the usual inflationary mechanism. We explicitly 
evaluate all the constraints on the model, under both the assumptions of prompt and delayed re- 
heating, and outline the viable parameter space. 
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I. INTRODUCTION 

There has recently been renewed interest in an al- 
ternative inflationary mechanism for generating an ap- 
proximately scale-invariant adiabatic density perturba- 
tion spectrum to the usual one. Rather than immedi- 
ately generating a curvature perturbation via perturba- 
tions in the inflaton field [1], instead the mechanism re- 
lies on isocurvature perturbations in another scalar field 
whose energy density is subdominant during inflation. 
After inflation ends this second scalar field comes to con- 
tribute significantly to the energy density, at which point 
the isocurvature perturbation converts to adiabatic even 
on super horizon scales. Subsequent complete decay of 
this second field guarantees purely adiabatic perturba- 
tions, though variants on this scenario can leave a resid- 
ual isocurvature component too. This mechanism of con- 
version from isocurvature to adiabatic perturbations was 
first discussed long ago by MoUerach [2] , was briefly men- 
tioned by Linde and Mukhanov [3] in a paper primarily 
focussing on scenarios for nongaussian isocurvature per- 
turbations, and more recently received renewed attention 
in Refs. [4, 5, 6]. Amongst these, Lyth and Wands [5] con- 
sidered the scenario in the broadest context, and named 
the second scalar field the curvaton. 

While Lyth and Wands described the development of 
the perturbations in considerable detail, they sought to 
keep their discussion as model-independent as possible 
and did not discuss a specific inflationary scenario. In 
this paper we give a specific realization of the curva- 
ton scenario and to evaluate all the constraints on model 
building that need to be satisfied for a successful scenario. 



II. THE SIMPLEST CURVATON MODEL 

While the general curvaton scenario allows perturba- 
tion generation featuring possible isocurvature compo- 
nents and possible nongaussianity, our aim here is to 
construct a simple curvaton model which creates a nearly 
scale- invariant spectrum of gaussian and purely adiabatic 
perturbations. The simplest possible curvaton model fea- 
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tures two massive non-interacting scalar fields, giving a 
potential 

y(0,a) = iAfV + ^™V, (1) 

where we indicate the inflaton by (j) and the curvaton by 
a. The scenario requires that the curvaton energy den- 
sity contributes negligibly during inflation (in particular 
the final stages). Since generically the curvaton will end 
up being the lighter of the two fields, it turns out that 
the curvaton must be close to its minimum, in order to 
prevent it driving a second period of inflation after the 
0-driven inflation is complete. Under these constraints, 
a will remain constant to a good approximation during 
the later stages of inflation, and we indicate this constant 
value by cr* . This initial condition is not fixed by the the- 
ory, but rather represents an additional free parameter to 
be fixed by observations.^ The curvaton's subdominance 
requires rn^al ^ M^0^ during inflation (where we will 
have (j) > mpi where wpi is the Planck mass). 

The early stages of inflation set the global mean of 
(7* in our observable Universe and arrange its classical 
homogeneity, and it then receives perturbations via the 
usual quantum mechanism, with the typical perturbation 
accrued in a Hubble time being 5a ~ H/2t:. In order for 
the eventual curvature perturbation to be gaussian, this 
perturbation must be small compared to the mean value 
of the field, and so we require al ^ H'^/Att'^. For these 
quantum perturbations to be the dominant influence on 
the curvaton, it must be effectively massless during infla- 
tion which requires <C — AirM'^cfP' /irn^y 

Once inflation is over and the inflaton energy density 
converted to radiation, the a field will continue to re- 
main constant while its mass is negligible compared to 
the Hubble parameter. Once ~ , it will begin to 
oscillate about the minimum of its potential, its energy 



Under this condition the oscillation mechanism of Ref. [7] is 
highly suppressed, and no mixture of correlated adiabatic and 
isocurvature perturbations will be relevant at the end of infla- 
tion. 
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density decaying at an average rate of p„ oc (the 
Universe wiU have to stiU be radiation dominated by this 
stage, as otherwise the domination of cr will initiate a 
new period of inflation). 

The final stage is the decay of the curvaton, which 
in this paper we will assume is a complete decay into 
conventional matter which thermalizes with the existing 
radiation. The decay occurs on a timescale r^, which 
is a further free parameter of the scenario. The most 
conservative constraint on the decay rate is that conven- 
tional radiation domination had better be in place by 
nucleosynthesis, though baryogenesis scenarios are likely 
to require a much earlier decay. Decay will happen when 



H, so this requires > i?n 



10-^"mpi. The 



time until decay sets the magnitude of the curvature per- 
turbation generated, because it determines what frac- 
tion of the mean energy density comes from the curva- 
ton when it decays. We require to match the COBE 
normalization; this sets a minimum requirement on the 
size of the curvaton fluctuations because they must be 
large enough to generate the required perturbations in 
the limit where the curvaton field is completely domi- 
nant when it decays. The inflaton will also generate a 
curvature perturbation at some level, and while a mixed 
perturbation scenario is permitted we will only consider 
here the case where the inflaton-generated perturbation 
is negligible, which requires M <C 10~^mpi. 



III. MODEL CONSTRAINTS 
A. The case of prompt reheating 

In this subsection, we shall assume that after infla- 
tion ends reheating occurs promptly, with the inflaton 
decaying into radiation. In that case, no further param- 
eters arc necessary to specify the scenario; wc have four 
parameters which are the two masses m and M, the cur- 
vaton value during inflation cr*, and the curvaton decay 
constant 

In a quadratic potential, inflation ends by violation of 
slow- roll at 0end — wpi/V^, corresponding to a Hub- 
ble parameter iJ^^d ~ During the subsequent 
radiation-dominated era H"^ = H^nd'^tnd/ where a is 
the scale factor. 

The next event to take place is for the curvaton to be- 
come eff'ectively massive, = H'^. This happens when 
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In order to prevent a period of curvaton-driven inflation, 
the Universe must still be radiation dominated at that 
point, which implies a significant constraint 



0-* < ^ TOpi • 



(3) 



This is a substantial restriction amongst all the possible 
values that a might have taken (most of which would 



result in a long epoch of fi-driven inflation after the (p 
field has reached its minimum). 

The most important constraint on the parameters is 
the requirement of reproducing the observed perturba- 
tion amplitude. Denoting the ratio of the curvaton en- 
ergy density to that of radiation by r = Pa/prad, in the 
limit where r < 1 (i.e. the curvaton decays during ra- 
diation domination), Lyth and Wands [5] demonstrated 
that the spectrum of the Bardeen parameter P^, whose 
observed value is about 2 x 10~^, is given by 



decay -n* 



where 



Hi ~ 
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(4) 



(5) 



is the Hubble parameter when observable perturbations 
were generated, around 50 e-foldings before the end of 
inflation, and 
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Continuing to presume that the Universe is still radia- 
tion dominated at decay, this will be at = H'^ = 

^endaend/«^ giving 



and hence 



rdecay- 3 



(for rdecay < 1) ■ 



(7) 



(8) 
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We find that achieving the correct perturbation ampli- 
tude, in combination with other constraints, excludes all 
the regions of parameter space where the curvaton decays 
while still effectively massless. 

In the opposite regime, where rdecay exceeds one and 
so decay occurs after curvaton domination, this formula 
no longer holds and instead the perturbation produced 
becomes independent oiV„, being [5] 



1 m 



9 ttV^ ~ 4a2 



M2 

- TZa (fo"^ '■decay > 1) • 



(9) 



The two expressions agree at the transition rdecay ~ 1- 
Note that if this last expression is normalized as quoted 
above, the gaussianity condition a1 » H^/Att^ ~ is 
automatically satisfied. 

With four parameters to vary and only one equality, 
Eq. (8) or (9), imposed upon them, we expect quite a bit 
of freedom in choosing suitable parameters. However the 
set of inequalities the parameters must satisfy is a large 
one, and it turns out that viable parameter space is quite 
restricted. For definiteness, we set thresholds that the in- 
fiaton perturbation must be no more than ten percent of 
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the curvaton perturbation (i.e. M < 3 x 10~^mpi) and 
that the gaussianity condition a J ^ HI/A-k"^ on observ- 
able scales be satisfied by an order of magnitude. 

We proceed by considering triplets of values 
(m, M, (7*). For such a triplet, we fix T„ using the power 
spectrum normalization. In some areas of parameter 
space this cannot be achieved, cither because the cur- 
vaton perturbations are too small even if the curvaton 
energy density becomes dominant, or because one would 
violate ^ 10~^°mpi as required by nucleosynthesis. 
Otherwise, we then test whether the many other require- 
ments to build a successful model are achieved, namely 
we must guarantee that the inflaton dominates during 
inflation, that the curvature perturbation from the infla- 
ton is negligible, that the curvaton is effectively massless 
during inflation and that the curvature perturbations re- 
sulting from the curvaton are gaussian, and that the cur- 
vaton energy density is still subdominant when it begins 
oscillating. This set of constraints slices off regions of 
the parameter space, leaving the region in which viable 
models can be constructed. 

Fig. 1 shows the allowed regions for two choices of cr*, 
with the main constraints plotted. The perturbation am- 
plitude constraint has two branches; the vertical part in- 
dicates that the models cannot reach the required pertur- 
bation amplitude even once the curvaton becomes fully 
dominant,^ while the lower part of the curve indicates 



^ On the vertical line itself lie models where the curvaton can dom- 
inate at decay; these are best analyzed separately which we do 
in the following paragraphs. 



that the perturbations have not grown sufficiently by nu- 
cleosynthesis. For low (7*, the nongaussianity constraint 
sweeps leftwards across the allowed region and there are 
no viable models once cr* < 10~"'^''mpi. For high cr* it 
becomes impossible to generate sufficient perturbations, 
cutting off parameter space above cr* ~ 3 x 10~'^mpi. 
There is however a significant parameter space of viable 
models satisfying all our requirements. 

For the special case of the curvaton decaying when its 
energy density is dominating over radiation, an analytical 
analysis of the parameter space is possible. In this case 
Fcr no longer enters into the relevant expression for T'^, 
Eq. (9). Thus, using the power spectrum normalization 
to fix the value of cr* in terms of M as cr* = 1.2 x 10^ M^, 
the relevant constraints reduce to 

M < 3 X 10"'^mpi, m^«;— (10) 



and 



ifn^ei ^ 10-4°mpi < F,, < m (^^y . (11) 



The first two conditions require that the inflaton- 
generated perturbations be subdominant and that the 
curvaton be massless during inflation respectively, and 
are the same as before. If they are satisfied, the remain- 
ing constraints automatically follow. The condition on 
Fo- derives from requiring rdecay > 1; where in this case 
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as a consequence of the a decay after curvaton domina- 
tion. 
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FIG. 2: The parameter space in the case of the curvaton 
dominating before decay. Models within the triangular region 
are viable; they have cr» = 1.2 x 10*^ M^, and may take on 
any value in the range given by Eq. (11). 



The formulae Eqs. (10) and (11) give an allowed re- 
gion in the m-M plane shown in Fig. 2. At each point 
within the allowed region there is a range of permitted 
r„ indicated by Eq. (11) which gives the correct density 
perturbation normalization. These regions correspond to 
the vertical part of the perturbation amplitude constraint 
in Fig. 1 for the corresponding cr*. 



B. The case of prolonged reheating 

Reheating is not expected to be instantaneous, and in 
this subsection we generalize the previous results to allow 
for a delay before reheating is complete. However there 
are no qualitative differences to the scenario and so we 
will keep the discussion fairly brief. 

In the case of prolonged reheating, after inflation there 
is a significant period during which the inflaton oscillates 
coherently at the bottom of its potential. Its ultimate 
decay products will be considered much lighter than (f) 
itself thus constituting the radiation, and we assume that 

decays into radiation with a rate T^. We make the 
assumption that there are no significant decays of the 
inflaton into the curvaton field. ^ The new parameter Fa, 



^ Were there such decays, they would typically generate an extra, 
nearly homogeneous, component of the curvaton energy density, 
which might affect the ability of the curvaton to generate suffi- 
ciently large perturbations. 



enters to modify the previous constraints. 

To prevent a new period of inflation driven by the cur- 
vaton, its energy density must still be subdominant dur- 
ing reheating. Following some standard approximations 
[8] , we can say that the inflaton (f) starts oscillating at the 
end of inflation (when i?end = and it behaves as 

P0 DC o~'^ until H fa T4,, when it decreases exponentially 
(/90 oc e~'"**) producing most of the radiation. When 
H fa reheating is completed. Before that point some 
radiation will be produced, but it is subdominant and 
Prad scales as a~^/^. After reheating the universe will be 
radiation dominated. 

Concerning the curvaton we can consider two situ- 
ations, one where a starts oscillating after reheating, 
and the other where it begins oscillating during reheat- 
ing. Just for illustrative purposes we assume that a fi- 
nally decays with a rate T„ after reheating when it is 
still subdominant with respect to the produced radia- 
tion. Then is given by Eq. (4) with rdccay < 1- The 
constraints defining the curvaton model during inflation, 
as described in Sec. II, still hold as does the require- 
ment Va- > i?nuci — 10~*°mpi. When cr starts oscillating 
it should be a subdominant component in order to pre- 
vent a period of curvaton-driven inflation, and it can be 
checked that this constraint remains the same as Eq. (3). 

In the case where the oscillations of the curvaton start 
after reheating > m and 
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fT' 



(13) 



where aroii is the scale factor at the end of reheating. 

On the other hand the produced radiation has an en- 
ergy density p\'^^ = SmpiF^/STt. As a consequence the 
expression for rdccay is exactly the same of Eq. (7), and 
it is independent of F^ 
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(14) 



Note, however, that if F^ > m and the constraint in 
Eq. (3) is satisfled, then 
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m 



(15) 
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This amounts to saying that throughout reheating the 
energy density of the curvaton has to be much smaller 
than the radiation energy density at the end of reheating. 
The strong constraint in Eq. (15) implies that the ratio 
r at the moment of the curvaton decay can indeed be 
much smaller than in the case of prompt reheating. In 
fact this case is recovered when F^ = H^^^ = M^/3, but 
for a prolonged period of reheating F^ < H'^^^^- Thus a 
similar analysis of the parameter space could be done as 
in the previous subsection, taking into account now that 
m2 < F2 < M2/3. 
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The curvaton begins oscillating during reheating when 

3 



^end 
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In the approximation ~ P^"'^(aend/areh)^, one finds 



^decay — 



in a 




(17) 



Since < m, a comparison with Eq. (14) shows that 
''decay Can be even smaller than in the case when a starts 
oscillating after reheating: once cr* , and are fixed, 
the curvaton energy density starts decreasing as at 
an earlier time. 



IV. CONCLUSIONS 

The curvaton model is an interesting new proposal for 
generating the approximately scale-invariant curvature 
perturbations which presently give the best match to ob- 
servational data. While nongaussianity would not nec- 
essarily be an observational disaster, we have chosen to 
restrict our attentions to gaussian models. We have con- 
structed the simplest possible realization of this idea, and 
demonstrated how various requirements close off regions 
of parameter space while leaving a substantial area of 
viable models. Even this simplest model, with prompt 
reheating, features four parameters, and so the predicted 
perturbations in different regions of parameter space are 
highly degenerate. 

While three of these parameters (to, M and F^) are pa- 
rameters of the underlying theory, the fourth, cr* , refers 
to the initial conditions for our patch of the Universe. 
The required values of ct* for a successful curvaton model 
have magnitude less than about 10~^topi, which repre- 
sents only a small region of the plausible values for ct* 
that might exist during inflation, as its potential has too 
low a magnitude to influence the dynamics. Nevertheless, 
in a typical chaotic inflation scenario one expects regions 
which do satisfy this criterion by chance. If one wishes, 
there are also opportunities to introduce anthropic prin- 
ciple considerations; for fixed values of to and M, regions 
with large a* would typically lead to the inflaton evolv- 
ing to the bottom of its potential followed by a period 



of slow-roll inflation driven by the curvaton, which given 
the small value of to would generate a much lower level 
of curvature perturbations than in the curvaton domain 
and hence not give rise to structure in the Universe. 

As far as the density perturbations are concerned, the 
predictions of these models are indistinguishable from 
slow-roll inflation models arranged to give the same spec- 
tral index. However, the curvaton model has the fea- 
ture that the gravitational wave amplitude is predicted 
to be low, because the Hubble rate during inflation is less 
than in an equivalent slow-roll model. In particular, the 
scalar and tensor perturbations typically will not obey 
the usual consistency relation (see e.g. Ref. [1]). In fact, 
in a natural curvaton model where a is extremely light 
(to^ <C H^) and subdominant with respect to the infla- 
ton, the scalar and tensor indices are predicted to be the 
same: Ug — 1 = rit + 2m^/3H^ ~ [5]. Equal spectral 
indices is a prediction of the power-law class of conven- 
tional inflation models, but they predict a high amplitude 
of gravitational waves, contrary to the curvaton model. 
Unfortunately however the low amplitude of tensors in 
the curvaton model makes it difficult to detect the ten- 
sors at all, and one cannot expect a useful measure of 
their spectral index. 

While one can certainly design slow-roll inflation mod- 
els with any spectral index for the density perturbations 
and negligible gravitational waves, so that the curvaton 
model predictions cannot be viewed as distinct, a de- 
tection of a signiflcant amplitude of gravitational waves 
would be sufficient to rule out these curvaton models. It 
is however interesting to note that present observations of 
large-scale structure and cosmic microwave background 
anisotropics tend to prefer a slight red tilt and a low con- 
tribution of tensor modes [9], as predicted by models of 
the type we have discussed. 
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